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We predict that the flow of unpolarized current in electron-doped GaAs and InP at room temper- 
ature is unstable at high electric fields to the dynamic formation of spin-polarized current pulses. 
Spin-polarized current is spontaneously generated because the conductivity of a spin-polarized elec- 
tron gas differs from that of an unpolarized electron gas, even in the absence of spin-orbit interaction. 
Magnetic fields are not required for the generation of these spin-polarized current pulses, although 
they can help align the polarization of sequential pulses along the same axis. 



(N 



c3 



-a 

o 
o 



> 

m 
o 
o 

-I— > 
c3 



i 

C 

O 

o 



X 



Spin-based semiconductor electronics ( "spintronics" ) 
requires the generation of spin polarization in nonmag- 
netic semiconductors 0, 0- Approaches through the in- 
jection of spin-polarized electrons or holes from magnetic 
materials yl Iftt la l a ML a nd spin- filteredjl or spm- 
selective |l(l 111! Il2t llfl. Il4j currents using the spin-orbit 
interaction have been considered. These all rely on an 
energetic difference between spin-up and spin-down car- 
riers in the electronic structure of some part of a circuit. 
Here we identify an entirely different approach. 

We show that initially unpolarized electron current 
flow in semiconductors can be unstable towards the spon- 
taneous formation of spin-polarized current pulses even 
without an applied magnetic field. The mechanism we 
propose is related to the charge Gunn effect [lf5j. in which 
homogeneous charge current flow is unstable to the spon- 
taneous formation of high-electric-field domains, inho- 
mogeneous charge distributions, and current pulses. We 
predict the "spin Gunn effect" can be seen at room tem- 
perature and pressure in GaAs and InP, and possiblv|l(ij 
GaN. A dependence of the electron drift velocity on spin 
polarization — originating from the Pauli exclusion prin- 
ciple — drives the spin Gunn effect. A room-temperature 
source of spin-polarized electrons from a nonmagnetic 
semiconductor would significantly advance the range of 
realizable spintronic devices. The high-frequency oscil- 
latory nature of these pulses also suggests new possibil- 
ities, such as may emerge from matching the oscillation 
frequency of the spin-polarized source to the precession 
frequency of spins in another material. 

The instability of the spin Gunn effect originates from 
an electron velocity (v = [iE, where \i is the mobility and 
E the electric field) that (1) depends on the local spin 
polarization of the electrons, and (2) differs for spin-up 
and spin-down electrons. The spin-dependent mobility 
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Most semiconductors have a ^ 0. The common scat- 
tering processes for electrons are independent of the spin 
of the electron, but do depend on its energy |l7j. As 
shown in Fig. ^ due to the Pauli exclusion principle the 
chemical potentials of spin-up and spin-down electrons in 
a spin-polarized electron gas differ from each other. As 
a result, for a degenerate electron gas the energy (and 
velocity) distribution of mobile spin-up electrons differs 
from that of spin-down electrons. Therefore the mobility 
of electrons for P ^ will differ from that when P = 0. 
For scattering of conduction electrons from ionized impu- 
rities, or acoustic phonons via piezoelectric coupling, or 
longitudinal phonons via Frohlich coupling (LO-phonon 
scattering), the low-energy electrons are scattered more 
than the high-energy electrons [l7j. as shown in Fig. ^ 
In the spin-polarized electron gas the increased energy 
of the spin-up electrons leads to longer scattering times 
and thus a higher mobility. The decreased energy of the 
spin-down electrons leads to shorter scattering times and 
a lower mobility. This situation corresponds to a > 
in Eq. If the dominant scattering process involves 
acoustic phonons coupling via the deformation potential 
(DP-phonon scattering), then the spin-up electrons scat- 
ter more than the spin-down electrons [17} , as also shown 
in Fig. ^ and a < 0. The only common scattering pro- 
cess that produces a = is neutral-impurity scattering. 
Thus the mobility of electrons in a spin-polarized semi- 
conductor almost always depends strongly on the spin 
orientation, even when there is no explicit spin depen- 
dence in the scattering processes themselves. At room 
temperature, and the elevated temperatures involved in 
the Gunn effect, LO-phonon scattering is likely to dom- 
inate, and we expect a > 0. Calculations for GaAs and 
InP for P -> of 
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is the spin polarization in the semiconductor, and a de- 
pends on the dominant spin- conserving electron scatter- 
ing mechanism. 



as a function of density at 300K and 500K, and as a 
function of temperature for n = I0 18 cm -3 are shown in 
Fig. |21 The spin-polarization-dependence of the mobil- 
ity increases with increasing density, and decreases with 
temperature, although \a\ > 0.1 for all of these scattering 
processes for temperatures below 500K at n = 10 18 cm~ 3 . 
We now consider the effect of a ^ on the spin polar- 
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FIG. 1: The electron sea for spin-up and spin-down electrons 
in an electron-doped semiconductor, (a). The sea is not spin 
polarized, and the mobility of spin-up and spin-down electrons 
is identical, (b). The sea is spin polarized, and the mobility 
of spin-up electrons is either smaller (as in the case of DP- 
phonon scattering) or larger (as in the case of LO-phonon, 
ionized-impurity, or piezoelectric-coupled-phonon scattering) 
than the mobility of spin-down electrons. In the first case 
a < and in the second a > 0. 
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FIG. 2: (ab). a as a function of conduction electron den- 
sity for GaAs and InP at 300K (dashed line) and 500K (solid 
line), (cd). a as a function of temperature for GaAs and 
InP for n = 10 18 cm -3 Green curves are LO-phonon scatter- 
ing, blue curves are ionized impurity scattering, red curves 
are piezoelectric-phonon scattering, and black curves are DP- 
phonon scattering. 



ization of electrons in an initially unpolarized Gunn do- 
main. The Gunn domain itself moves in a self-sustaining 
way — an inhomogeneous electric field moves the charge, 
and the shifting space-charge region moves the electric 
field — and the net result is the collective motion of the 
electric field shape and the space-charge region through 
the sample. In the presence of these electric fields a 
small difference in mobility between spin-up and spin- 
down electrons means fewer electrons of one spin species 
will flow in response to the field. This enhances the dif- 
ference in mobility between spin-up and spin-down elec- 
trons, which leads to further differences in the accumula- 
tion or depletion of the two spin species. A positive feed- 
back effect therefore drives the increase in the spin po- 
larization. The electric field is determined by the charge 
density through the Poisson equation, and does not de- 
pend on the spin polarization, so the domain continues 
to move as before but now with an inhomogeneous spin 
polarization. We require only a small initial imbalance 
between spin-up and spin-down electrons to begin the ef- 
fect. This can be nucleated thermally as an initial inho- 
mogeneous spin polarization with a random orientation. 
It can also be nucleated with a small applied magnetic 
field, which will also serve to align the spin polarizations 
of all the domains. 

As the electron momentum scattering time is very 
short (< 100 fs) we can calculate the time dependence 
of the spin polarization using spin-resolved drift-diffusion 



equations^], 
~~dt 

at 



2Ti dx 



D 



2Ti 



dx 



D 



<9 2 rij 
dx 2 
d 2 n l 
dx 2 



(4) 
(5) 



where T\ is the spin relaxation time of the electrons. Us- 
ing Eq. QJ, keeping only terms up to first order in P, 
neglecting pure diffusion effects and shifting to a moving 
reference frame with the domain, x' = x — Vd om t , where 
Vdom is the velocity of the domain, 

dP 
~dt 
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The form of this equation yields spin amplification if 
7 > . The dP/dx' term describes the flow of in- 
homogeneous spin polarization from one region of the 
domain to another; similar term s ap pear in the charge 
Gunn effect itself and are ignored[ljj — they change the 
quantitative results somewhat but not the qualitative be- 
havior of the Gunn instability. In our treatment here we 
will neglect this term, which we calculate to be over an 
order of magnitude smaller than jP. The central quan- 
tity is the spin amplification rate 7. 

The first observation about Eq. J7J is that 7 is propor- 
tional to the spatial variation of the drift current. The 
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drift current varies spatially due to the inhomogeneous 
electric fields and densities associated with the Gunn do- 
main, and its derivative can be positive under certain 
conditions for a > (e.g. downstream from the domain 
center) and under others for a < (e.g. upstream). 
From the continuity equations Eqs. I|4I5I) . dn s /dt — Vj s , 
so differing spin-up and spin-down currents (from differ- 
ent mobilities) produce different spin-up and spin-down 
density accumulation. Far away from the domain the 
local charge density is the equilibrium charge density, 
and 7^0. To make 7 as large as possible we should 
consider cases involving large charge imbalances (within 
the breakdown fields of the material). Ordinary dop- 
ing levels for the charge Gunn effect are in the range of 
10 14 - 10 16 cm" 3 . Doping levels of 10 18 cm" 3 lead to 
sharper domains and less microwave power in the funda- 
mental mode of the oscillation, so they are not commonly 
used. For a spin Gunn effect, however, a sharp domain 
will lead to larger spin polarizations. 

We now quantitatively calculate the value of 7 for a 
doping level of 10 18 cm" 3 at a 300K lattice temperature 
in GaAs and InP. We ignore higher-order effects and eval- 
uate — (a/ n)d[nE fi(n / 2)] / dx' for P = 0. Hence 



7 = —afj,E 
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The first term above comes from the derivative of the 
electric field, which is related to the local charge density 
through the Poisson equation. For 7 to be positive, if 
a < the charge density must be larger than equilibrium 
(n > n ), whereas if a > the charge density must 
be depleted (n < n a ). Evaluation of 7 thus requires 
calculation of n(x') and E(x') for the Gunn domain. 

Our calculation of n(x') and E(x') follows that of 
Szc 20] for a mature, steady-state domain. We assume 
the electrons in the lower valley and the upper valleys 
have the same temperature 0, |2^ , and that the diffu- 
sion constant D is independent of the electric field, the 
electron density and spin. The time-dependent electron 
current is equal to the displacement current, 
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For a high-field domain propagating without a change of 
shape, in the moving frame (x' = x — Vdomt) 
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The electron drift velocity outside of the domain, vr — 
J/{eno)- The solution of equation (|TU|l is: 
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A self-consistent solution to Eq. (|11|) produces the re- 
quired relationship between the density and the electric 
field E(n). The spatial and temporal behavior of E can 
then be found by evaluating Poisson's equation, 
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where x' is the location of the peak electric field in the 
moving frame of reference. The local electron tempera- 
ture T e is determined by the local electric field from 

k B T e = k B Ti + (2/3)qT E [iE 2 [l + Rexp(~AE/k B T e )}- 1 

(13) 

where T; is the lattice temperature, te is the energy re- 
laxation time (~ 10" 12 s), AE is the energy separation 
between the T and L valleys, and R is the ratio of the 
density of states in the L valleys to the T valleys. We 
have used values for these quantities from Sze[2fJ. Solu- 
tions of Eqs. Hllfl - i|12|) f° r the density and electric field 
as a function of position are shown in Fig. for GaAs 
and InP for a 300K Tj. Fig.^Jab) show the density and 
Fig. EI cd) show the electric field. 

The amplification factor 7 is shown as a function of 
position in Fig. OHef ) for the two materials if the domi- 
nant scattering process is DP-phonon scattering, and in 
Fig. Elgh) for LO-phonon scattering. As the electron 
temperature varies with position according to Eq. Ijl3|l . 
a must be evaluated as a function of temperature and 
density to determine 7 properly. The electric field, and 
thus the spatial variation of the drift current, is largest 
at the center of the domain, however the electron tem- 
perature is also greatest there. These competing effects 
produce a maximum 7 which is near the domain but not 
at the peak of the space charge. For GaAs there are re- 
gions of amplification for spin lifetimes as short as 3 ps, 
whereas for InP the amplification is present for lifetimes 
longer than 2 ps. The region of spin amplification differs 
for the two scattering mechanisms; for LO-phonon scat- 
tering the largest spin amplification is downstream of the 
domain, whereas for DP-phonon scattering the largest 
spin amplification is upstream of the domain. 

The significant increase in 7 with doping density leads 
us to propose Gunn diodes with n ~ 10 18 cm" 3 for the 
spin Gunn effect. For a lower doping density of 10 16 cm" 3 
the Ti of GaAs at 300K is 50 ps|23j and our calculated 
7 ~ 4 ns , so spin amplification is not expected to oc- 
cur. The elevated temperature near the Gunn regime, 
which ranges from 500K a distance 100 nm away from 
the center of the domain, to much higher values near the 
peak electric field, will further reduce the spin relaxation 
times. Scaling with T 3 from D'yakonov-Perel' preces- 
sional relaxation |23j suggests spin relaxation times four 
times shorter at 500K, of the order of 12 ps. The reduced 
mobility at higher temperatures will lengthen these esti- 
mated times, and our calculated spin amplification rates 
of j(P — 0) > 0.4 ps" 1 are five times larger than this 
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FIG. 3: The spatial dependence of the electron density n 
[(a) for GaAs, (b) for InP] and the electric field E [(c) for 
GaAs, (d) for InP] in the frame moving with the Gunn domain 
(x' = is the center of the domain). The spin amplification 
rate 7 for DP-phonon scattering [(e) for GaAs, (f) for InP] 
and for LO-phonon scattering [(g) for GaAs, (h) for InP]. The 
three curves correspond to different drift velocities far from 
domain (different applied voltages to the Gunn diode): black 
is 5.7 x 10 6 cm/s for GaAs, 2.6 x 10 7 cm/s for InP, red is 
5.0 x 10 6 cm/s for GaAs, 2.45 x 10 7 cm/s for InP, and blue is 
4. Ox 10 6 cm/s for GaAs, 2.4x 10 7 cm/s for InP. The saturation 
polarization (corresponding to j(P) = T^ 1 ) is shown in (i) 
for GaAs and (j) for InP. 



rapid spin relaxation rate, providing confidence that spin 
amplification is likely for room-temperature devices. The 
enhanced electron-electron scattering present in degen- 
erate systems at these high temperatures is expected to 
further increase the spin relaxation times (without af- 
fecting the mobility) |2J] . Shown in Fig.^ij) are steady- 
state (saturation) values of the spin polarization (P sa t) 
estimated for LO-phonon scattering determined by set- 
ting T\ — 10 ps, using a(P ^ 0) in Eq. JSJ, and solving 
7(P) = Tf 1 . As a(P = 1) = 0, a changes 100% between 
P = and P = 1. In contrast, the total n, total /i, and E 
change ~ 20%. Thus we neglect nonlinear effects on the 
total n, total [i, and E in our calculation of P sa t- When 
7(P = 0) < T± there is no amplification and P ea t = 0. 
The largest P sat exceeds 80% for both GaAs and InP, 
and should be directly visible in a Faraday rotation mea- 
surement (e.g. Ref. 25]). Therefore the spin Gunn effect 
is a pulse of highly spin polarized electrons located just 



before or just after a charge current pulse. 

The existence of this spin Gunn effect is exceptionally 
robust to temperature and spin relaxation, suggesting 
a wide range of potential applications of devices based 
on the effect. We have only considered here spin effects 
amenable to nearly analytic analysis, such as occur in the 
propagation of a "mature domain" . Other known modes 
of operation of the charge Gunn effect, such as the lim- 
ited space-charge accumulation (LSA) mode, will provide 
additional spintronic functionality when analyzed in the 
context of the instability to non-zero spin polarization we 
have identified here. 
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